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Abstract 

We investigate two-photon double ionization of helium by intense (~ 10 15 W/cm 2 ) ultrashort 
(~ 300 as) soft X-ray pulses (E = 91.6 eV). The time-dependent two-electron Schrodinger equa- 
tion is solved using a coupled channel method. We show that for ultrashort pulses the angular 
distribution of ejected electrons depends on the pulse duration and provides novel insights into the 
role of electron correlations in the two-electron photoemission process. The angular distribution at 
energies near the "independent electron" peaks is close to dipolar while it acquires in the "valley" 
of correlated emission a significant quadrupolar component within a few hundred attoseconds. 

PACS numbers: 32.80.Rm, 32.80.Fb 
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I. INTRODUCTION 



Recent advances in the high-order harmonic generation (HHG) techniques have led to the 
development of soft X-ray sources that feature ultrashort pulses with pulse durations of a few 
hundred attoseconds (as) jl] and may reach intensities (> 10 14 W/cm 2 ) that are capable of 
inducing multiphoton processes. Extreme ultraviolet (XUV) pulses (photon energy 27.9 eV) 
with pulse duration of 950 as have been characterized with an autocorrelation technique 
Q]. Recently, the two-photon double- ionization and above-threshold ionization of helium 
were experimentally observed with the Ti:sapphire 27th harmonic pulses (photon energy 
41.8 eV) |3[. These experimental advances open up the opportunity to revisit the dynamics 
of double ionization of helium by XUV photons previously investigated only in the single- 
photon absorption and scattering regime using synchrotron radiation jj, |^. Simultaneous 
ejection of two electrons by a single photon allowed detailed tests of wavefunctions for 
the three-body Coulomb problem jy, Q, [J and the role of electron correlations in strongly 
inelastic processes accompanied with near-zero momentum transfer (photoabsorption) or 
sizable momentum transfer (Compton scattering) P. [lo|. 

Multi-photon, in particular two-photon, ionization of helium by XUV pulses has been 
studied theoretically by different groups. A considerable numerical effort has been made 
to solve the two-active electron time-dependent Schrodinger equation (TDSE) with various 
methods. The R-matrix Floquct theory was successful to describe the (27,2c) process of 
He [lj] in the photon energy range between where absorption of two photons are necessary 
for double ionization. The configuration interaction B-spline spectral method jl^ . Ik^ was 
applied to solve the TDSE for this problem. The products of two B-splines represent the 
radial part of the wavefunction which allows the inclusion of the electron-electron interaction 
to a high degree of accuracy. Colgan et al. jlj] developed a time- dependent coupled chan- 
nel method and studied the complete fragmentation of helium at 45 eV photon energy and 
presented fully differential cross sections. Recently Lambropoulos et al. |l5j] found a "knee" 
structure in the intensity dependence reminiscent of a similar knee shape for double ioniza- 
tion by strong IR pulses [16(. Photons above the double ionization threshold (o; XU v > 2.9 
a.u. or 79 eV) were considered by Parker et al. |l7j who performed the direct numerical 
integration of the two-electron TDSE with a mixed finite-difference/basis set approach on 
a lattice and studied double-ionization with 87 eV photon energy pulses with a laser peak 
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intensity around 10 16 W/cm 2 . They analyzed both sequential as well as non- sequential dou- 
ble ionization events by a varying number of absorbed photons for long pulses (r p > 2 fs). 
Most recently, Ishikawa and Midorikawa 18( investigated two-photon double ionization by 
ultrashort pulses with durations of r p « 150 to 450 as pertinent to HHG sources. They 
identified an "anomalous" component in the electron spectrum in between the two peaks 
associated with sequential double ionization and discussed its possible origin in terms of 
post-ionization energy exchange and core relaxation effects. 

In this paper, we theoretically investigate two-photon double ionization of helium by 
ultrashort attosecond pulses as a function of time by solving the TDSE with our coupled 
channel method which has been originally developed for heavy-ion helium collisions 
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ater implemented to describe laser-driven atomic processes and two-photon coherent 



control 22]. We consider experimentally realized high intensity laser pulses with 13.5 nm 

wavelength which are the 59th harmonic of a Ti:sapphire laser (wavelength 800 nm). 

The photon energy considered (91.6 eV) is larger than the double ionization threshold of He 

(79 eV). A single photon is thus sufficient to induce double ionization. This case has been 

studied in detail with weak- field synchrotron sources where multi-photon effects are absent. 

Because one photon can interact with one electron only, double ionization cannot occur 

without electron-electron interaction. The picture is that one electron is directly ionized by 

absorbing the photon, and the second electron leaves through electron correlation either in 

the initial or in the final state, or both. This has been discussed in terms of a shake-off and 

' I — I 

electron-electron scattering (often referred to as TS1 [24]). 

In contrast, for intense fields considered here, there is sufficient photon flux such that 
two photons can be absorbed, one by each electron, and the pair of electrons is ejected. 
Electron correlation is therefore not a prerequisite for double ionization to occur. At the 
same time, the ponderomotive energy of the XUV pulse E P = 8nl /4cu^ UY is so small that 
ionization by the rescattering of the first ionized and accelerated electron that causes "non- 
^uentiaT doub.e _ by strong IR Pul8 es 0H can be rul ed out. We d 18 en S8 the 
conceptual difficulties in applying notions of sequential and non-sequential double ionization 
to such short pulses. We show that the angular distribution provides detailed insights into 
the ionization process on the attosecond time scale. The role electron correlation plays in 
this process can be identified. 



3 



II. SCENARIO FOR TWO-PHOTON DOUBLE IONIZATION BY ATTOSECOND 
XUV PULSES: TIME SEQUENCE AND CORRELATION 



We consider a linearly polarized attosecond XUV pulse with a Gaussian envelope, 



F(t) =F exp 



-2 In 2 



T 2 
P 



cos(u; xuv t)e^ , (1) 



where t p is the full width at half-maximum (FWHM) of the pulse intensity. The center 
frequency w XU v corresponds to a Ti:sapphire 59th harmonic pulse with the energy of 91.6 
eV. Following Ref. 18] we will consider pulse durations t p = 150 as and 450 as corresponding 
to r p = 6.25 and t p = 18.7 a.u. The period of the XUV cycle is T = 2tt/ujxuv — 1-9 a.u. 
The XUV pulse (Eq. subtends only few cycles (3 to 10) and therefore closely resembles 
few-cycle optical or near-IR pulses. The significant Fourier broadening therefore precludes 
the appearance of spectrally sharp photoionization peaks. There is, however, a fundamental 
difference to optical pulses of the same intensity: even at an intensity of I = 10 15 W/cm 2 , 



the quiver amplitude of a free electron, ^8nl /c/uj\ nY ~ 0.01 a.u. is small on an atomic 
scale. Likewise, the ponderomotive energy Ep = 8nl /4cu^ nv ~ 0.0006 a.u. is negligibly 
small. Therefore, ionization takes place deep in the (multi) photon regime rather than in the 
tunnel ionization regime applicable to IR pulses of the same intensity. This difference has 
immediate consequences for the notion of "(non) sequential" ionization. While for tunnel 
ionization the time window At of an individual ionization "burst" can be uniquely identified 
near the field maxima with sub-cycle precision 27], the multi-photon ionization event is 
intrinsically delocalized in time over several cycles. Only then does the electron response 
to an electromagnetic pulse mimic that of photon absorption. In view of the fact that the 
entire XUV pulse duration r p subtends only a few cycles, it is obvious that the notion of 
sequentiality of ionization events loses its meaning in the present case. This is in sharp 
contrast to intense field ionization by optical fields. There, the first ionization by tunnel 
ionization under the influence of a quasi-classical electric field is well localized and separated 
in time from the collisional ionization of the second electron upon rescattering. The observed 
scaling with the pulse duration oc Tp , where N is the number of photons absorbed, should 
therefore not be taken as evidence of (non) sequentiality but a measure of the total energy 
absorbed from the radiation field during r p . The uncertainty in time when the absorption 
process takes place or time derealization of the multi-photon processes does not imply that 
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all time-differential information on the ionization process is averaged out, as will be shown 
below. 

The time characterizing the pulse duration should be compared to the time scale of the 
electronic motion. Using the approximate hydrogenic expression for the classical orbital 
period 

r = 27in 3 /Z 2 eff , (2) 

the orbital period ranges between To = 40 as (= 1.6 a.u.) for the "inner" electron of 
He + (ls)(Z e ff = 2) and, for the "outer" electron of He (Is 2 ) with a binding energy of 
24.6 eV (Zgfj = 1.3), r = 90 as (~ 4 a.u.). The cycle period T and the orbital period 
are comparable to each other, thus probing the electronic motion on the time scale on 
which the two interacting electrons of the helium ground state exchange energy, linear and 
angular momentum. Thus, double ionization by attosecond XUV pulses may probe electron 
correlations in both initial and final states. 

The role of correlation in double photoionization of helium is a well-established subject in 
the low-intensity or single-photon limit of XUV radiation going back to the pioneering paper 
by Byron and Joachain 28|. As the electron-photon interaction is a one-body operator, 
single-photon absorption can directly eject only one electron. Ejecting a second electron 
requires with necessity electron-electron interaction. The latter does not, however, imply 
correlation effects. Adhering here and in the following to the identification of correlation 
with those pieces of the interaction not included in a mean-field or independent particle (IP) 
model as embodied in the (single configuration) Hartree-Fock description j^, one-photon 
double ionization can proceed via mean-field contribution. Already the sudden change of 
the screening following the ejection of the first electron generates a finite probability of 
ejecting a second electron. This "shake-off" process accounts for about 40 % of the total 
doble ionization cross sections at high photon energies. Clearly, for a quantitatively accurate 
description, in particular over the entire range of photon energies from threshold to high 
energies js| correlation effects beyond the mean field in both the initial and final states are 
essential. 

For two-photon double ionization by XUV pulses with Uxuv > 2 a.u. dominance of 
independent particle (IP) ionization is expected since each reaction 

He + hwxuv He + (nZ) + e~, (3) 
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FIG. 1: Electron spectrum following two-photon absorption {uj = w X uv = 91.6 eV) in coincidence 
with He 2+ , schematically. The symmetric energy is Es = Wxuv — {Ii + ^)/2 with I\2 being the 
first and second ionization potentials, respectively. 



where n and I are the principal and angular momentum quantum numbers of He + , respec- 
tively, is energetically allowed for all n. The quantum numbers, n = 1 and I = are expected 
to dominate in Eqs. Thus, correlation effects appear to be unimportant for the two- 

photon process. It should be noted that the picture of a chain of reactions each satisfying 
energy conservation in the photoelectric effect separately, invoked by Eqs. (0 EJ) is only 
meaningful for r p — ► oo. In this limit, Eqs. (13141) implies an electron spectrum in coinci- 
dence with He ++ (displayed schematically in Fig. 1) with two Rydberg series symmetrically 
centered around the energy 



The single-photon double ionization spectrum well-known from synchrotron-studies appears 
as a continuum below 12.6 eV (cj X uv — h — h)- Its two-photon replica would set in above 
E 

— ^xuv = 91.6 eV. For ultrashort r p all discrete peaks get dramatically broadened and 
merge into a quasi-continuum. 

Apart from the broadening, the limit of short r p has further consequences when this time 
becomes comparable to the electronic correlation time tc in the helium ground state which 




(4) 



£s = ^xuv-(/i + / 2 )/2 = 52eV 



(5) 
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can be simply estimated from the characteristic time for collisional exchange of energy and 
angular momentum between two classical electrons. Alternatively, it can be estimated from 
the correlation energy Eq — E — Ejjp as Tc = ^-/Ec- In either cases, tq is of the order 10 
a.u. (or 200 as). XUV pulses with periods T of 2 a.u. and durations of 3 - 10 a.u. therefore 
can probe the correlation dynamics. 

<a) (b) Mi MU 



FIG. 2: Schematic interaction diagrams for two-photon absorption from He ground state, (a) 
Each electron absorbs one photon each; (b) one electron absorbs two photons. Dashed lines denote 
the electron-electron interactions. 

It is instructive to visualize the two-photon double ionization process diagrammatically 
(Fig. 2). The two photon lines each representing the one-body operator of photoabsorption 
end either at the same or at two different electrons resulting in two different diagrams. (The 
line representing the nucleus has been omitted for simplicity). It should be noted that a 
definite time ordering of the vertices of electron-photon interactions is neither implied non 
meaningful for ultrashort pulses in light of the discussion above. The dashed lines refer to 
electron-electron interaction in the initial and final states which lead to energy and angular 
momentum exchange. The latter is reflected in a configuration-interaction wavefunction in 
terms of admixtures of orbitals of different single-particle angular momenta, 

1*0 = E i**> +E4 y V>b i > ( 6 ) 



fsee also 



Typical orders of magnitude of admixture coefficients for the initial state are 
Eq. (fTK|) below) \a p /a s \ < 0.1, \ad/a s \ < 0.01 and those of higher I are exceedingly small. 
The admixture of non-s orbitals to the He ground state provides a unique signature of 
electron correlation as it would be absent in an IP or HF model. More precisely, / ^ 
configurations represent angular correlation while coefficients a^'^ may contain radial 



correlation. Correspondingly, angular momentum components in the final state reflect both 
the angular momentum transfer Al = ±1 by photoabsorption as well as the non-s admixtures 
due to electron correlations in the initial and final state. Their presence can be mapped out 
by the time-dependence of the angular distribution of ejected electrons. 



III. COMPUTATIONAL METHOD 



We have calculated the double ionization by two-photon absorption represented by the 
diagrams of Fig. 2 using our time-dependent coupled-channel method. The point to be 
noted is that while we discuss and interpret our results within to lowest-order perturbation 
theory (LOPT), the calculation is fully non-perturbative taking into account electron-photon 
and electron-electron interactions to all orders, albeit within a truncated basis. Briefly, we 
calculate ionization process of the helium atom in the laser pulse by solving the TDSE 
equation 

z— vj/( ri , r 2 , t) = (H He + V{t)) *(n, r 2 , t), (7) 
for the atomic Hamiltonian, 

P? P? 2 2 1 

H B . = % + % + | r , 8 

2 2 ri r 2 |r x — r 2 | 

and with the laser-electron interaction 

V(t) = Y1 F ( £ ) • r * ( 9 ) 

i=l,2 

in the length gauge and the dipole approximation. The laser pulses are linearly polarized 
along the 2-axis with the time dependence given by Eq. (JTJ). We expand \t / (r 1 ,r 2 ,t) in the 
basis of eigenfunctions of the time- independent Schrodinger equation 

#$,(r 1 ,r 2 ) = E i $ J (r 1 ,r 2 ), (10) 

to yield 

N 

*(r ls r 2 ,t) =£a i (t)$ i (r a ,r 2 )e- iE ' t , (11) 

where the a,j(t) are the time-dependent expansion coefficients and Ej are the eigenvalues 
of (Eq. [TDJ). Inserting (Eq. ITTjl into the TDSE (Eq. |7j) leads to the system of first-order 
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differential equations for the expansion coefficients 



N 



dt 



(12) 



Denoting the ground state by k — 1, we impose the initial condition 



1 



k 



1 



a k (t 



oo) 



(13) 







k^l. 



The asymptotic probabilities for transitions into final states k after the pulse has been turned 
off are given by 



The ionization probability can be retrieved from P k which includes discretized channels 
representing the continuum formed by the wave packets. The equations of coupled channels 
(Eq. IT2"j) are solved by a Runge-Kutta-Fehlberg integrator of order five with automatic time 
step adjustment. 

The eigenfunctions $j in (Eq. nT)|) are obtained by diagonalizing the Hamiltonian in a 
basis of orthogonal symmetrized two-particle functions 



In the following we restrict ourselves to singlet helium states only. The two-particle functions 
are made up of symmetrized single particle orbitals, g e i{r)Y[ m , where the radial functions g e i 
consist of radial Slater functions and radial regular Coulomb wave packets. We note that 
the coefficients b$ are related to the admixture coefficients discussed earlier following Eq. 
(JHJ). The wave packets form a discrete representation of the Coulomb continuum and can 
serve as building blocks of our finite basis [3 l^tj] . 

We include single-particle wavefunctions with < l\, li < 2 angular momenta and couple 
them to < L < 2 total angular momentum two-electron states. For the L = configura- 
tions we use ss+pp + dd angular correlated wavefunctions, for L = 1 we use sp+pd couplings 
and for L = 2 the sd + pp + dd configurations, respectively. Since already the contribution 
of d orbitals in the present case is found to be small, higher U can be safely neglected. The 
angular correlated contributions play an essential role to understanding the angular distri- 
bution of the ionized electrons. In order to determine the final electronic state population, 



P k = W{t^ +oo)| 2 . 



(14) 




(15) 
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the expectation value of the reduced one-electron density operator p = J2i=i 2 <K r — r i) * s 
calculated after the laser pulse, 

p(f) = (*(*-> 00) \p\V(t + 00)). (16) 

We employ the Feshbach projection method [2^ to separate the singly-ionized states from 
the doubly-ionized states. Accordingly, the one-electron polar angular distribution of ionized 
electrons in the double ionization channel is given by 

2lT OO 

Pdi(0) = J(^ DI \J2 6 ( r - r ^Di)r 2 drd^ 

l=1 ' 2 

2ir oo 

= IJJJ ^Di^e^-Y^dY^drdip. (17) 

ri 

where ^di represents the projection of \I/ onto the subspace of doubly ionized states. 



IV. RESULTS AND DISCUSSION 
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FIG. 3: Energy distribution of the ejected electrons in two-photon double-ionization of He. The 
peak intensity of the pulse is 10 15 W/cm 2 and the pulse duration is r p =450 as. The solid line 
represents our results and the dashed line represents the data for Ishikawa et al. |la |. The energy 
positions referred to in text are 37.2 eV (a), 52 eV (b), and 67 eV (c). 

Before analyzing the angular distribution from Eq. ([17)1 we briefly present results for 
the energy distribution for which a direct comparison with a recent calculation by Ishikawa 
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and Midorikawa [18| is possible. The single electron energy distribution integrated over the 
second electron for the pulse duration of 450 as (Fig. 3) features two prominent peaks which 
can be easily identified with help of Fig. 1 as the ionization spectra following the reactions 
Eqs. (j21and0J). The dominant yet strongly broadened peak at 67 eV (labeled c) is due to 
electrons ejected from the ground state of He with the first ionization potential of 1% = 24.6 
eV. In the second interaction, the electrons are ejected from the He + ion with an ionization 
potential of I? = 54.4 eV, yielding the peak at 91.6 — I2 = 37.2 eV (labeled a). From the 
higher members or the Rydberg series only n = 2 peaks are identifiable in Fig. 3 as local 
humps, one just below and one above the main peaks, respectively. Structures from n > 3 
are not visible since their contributions become exceedingly small. The cross section of the 
single-photon double ionization continuum below 12.6 eV (see Fig. 1) is by far too small to 
be visible on a linear scale. The peaks a and c (Fig. 3) have been previously referred to as 

n n 

sequential ionization [18|| or above-threshold ionization [17]. We will refer to this process 
as independent particle (IP) ionization to stress that electron correlation effects play no 
significant role in their occurrence. This is in striking contrast to the spectral feature in 
the "valley" (labeled b) (also referred to as anomalous component Q|) in which correlation 
effects are of crucial importance. We refer to this feature as "correlation induced" (CI) 
ionization. 



10 20 30 40 50 60 70 80 90 100 
Electron energy (eV) 

FIG. 4: Notation is the same as in Fig. 3, but for r p = 150 as. The marked energy positions are 
39 eV (a), 52 eV (b), and 67 eV (c). 
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The dependence of the CI ionization on the pulse duration t p is illustrated in Fig. 4 for 
an ultrashort pulse of 150 as. The valley is now quite shallow and a significant fraction 
of the ionization probability is contained in the "valley". This is, in the first place, an 
obvious consequence of the increased Fourier broadening in the ultrashort pulse limit. In 
the opposite limit r p ^> T and t p 3> r , the spectrum is expected to revert to the quasi- 
discrete line spectrum, schematically depicted in Fig. 1. Only in the long-pulse limit, the 
notions of time ordering or sequentiality of the process takes on physical meaning. Overall, 
our data agree with those of Ref. 3| remarkably well on an absolute scale with the largest 
discrepancies in the wings of the peak for the 10 cycle pulse (450 as). 



T p =450as 




T D =450as 




FIG. 5: The angular distribution (polar plot) of the ejected electrons for an XUV pulse with r p = 
450 as. Snap shots of lines of constant intensities are taken at times 50 to 450 in steps of 50 as 
(from inside going outward) after the pulse's rise to half maximum, for energies 37 eV (left) and 
52 eV (right). The unit circle indicates intensities of 7 x 10 -6 eV -1 for (a) and 6 x 10~ 7 eV -1 for 
(b). The arrows show the polarization axis. 

The identification of the valley region near E s with correlated ionization is, in the first 
instance, taken over from one-photon double ionization by synchrotron radiation where 
the region of symmetric energy sharing of the available photon energy is dominated by 



correlation effects 



24j. The extension of this identification to two-photon absorption can 



be quantitatively justified by the properties of the angular distribution, as shown below. 
Were the valley simply the result of the Fourier broadening of two IP peaks, the distribution 
Pdi{6) at the energies near (Es) should closely resemble those of the spectral regions (a) or 
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(c). That this is not at all the case is illustrated by the polar plot (Fig. 5) of the angular 
distribution near the IP ionization peak (a) and the CI ionization valley (b). The different 
contour line indicates the time evolution of the angular distribution in increments of 50 as 
for the 450 as pulse. While the IP peak retains the emission pattern of a Hertz dipole during 
the entire pulse duration, the CI electron distribution takes on a pronounced non-dipolar, 
i.e. quadrupolar, pattern after about 150 as. The onset of a non-dipolar distribution on this 
timescale can be also observed for the ultrashort pulse of r p = 150 as (Fig. 6) indicating that 
the sharp differences in the angular distribution between the IP peak and the CI component 
is also present when the valley is very shallow. For the ultrashort pulse a slight peak shift 
from 37 eV to about 39 eV is found in agreement with Ref. 



FIG. 6: Polar plots of the angular distribution of the ejected electrons after a pulse with r p = 150 
as. The distribution (a) is taken at 39 eV and (b) at 52 eV. The unit circle indicates an intensity 
of 3 x 10~ 7 eV -1 . The arrows show the polarization axis. 

The anisotropy of the angular distribution can be characterized by the multipole expan- 
sion 



where o"o is the integral cross section, P 2 ,4 ar e the Legendre polynomials, and (3 and 7 are 
the second order (k = 2) and fourth-order (k = 4) anisotropy parameters, respectively. 
Note that a "dipolar" emission pattern has k = 2, i.e., it represents "alignment", while 
the "quadrupolar" pattern is of rank k — 4 and should be more correctly referred to as 
"hexadecapole" . Individually, the range of the multipole parameters are — 1 < /3 < 2, 





(18) 
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and —1 < 7 < 7/3, the highest order of anisotropy k = 4 is consistent with two-photon 
absorption, 

k<2N. (19) 

Higher anisotropy coefficients beyond k = 4 are not detectable. By projecting the numer- 
ically calculated angular distributions to Eq. (fT8|) . we obtain the (3 and the 7 parameters 
listed in Table 1. Near the IP ionization peaks, (37 and 39 eV for 450 and 150 as pulses, 
respectively) (3 is at least one order of magnitude larger than 7. The /3 values are very close 
to their maximum value of 2. 

At 52 eV, near the CI valley, j3 and 7 become comparable, giving rise to a strong mixing 
of dipole and quadrupole terms in the angular distributions. 

TABLE I: The multipole expansion parameters (3 and 7 Eq. ()18j) for the pulse durations of 450 
and 150 as, at two energies each, corresponding to IP and CI ionization (see text). 

t p Type Energy (3 7 
450 as IP 37 eV 1.94 -0.08 

CI 52 eV 0.40 0.58 
150 as IP 39 eV 1.87 -0.17 

CI 52 eV 0.51 0.35 



It should be noted that the present deviation from a strictly aligned {k = 2) pattern is due 
;o multiphoton effects and not due to retardation effects beyond the dipole approximation 
29^ . Figs. 5 and 6 clearly show that two-photon IP ionization features a near Hertz dipole 
distribution while CI ionization possess a significant k = 4 admixture. It is now instructive 
to relate the origin of the quadrupole component to correlations. A non-vanishing 7 requires 
a final state in the continuum with Lf = 2 since k = 2Lf. The latter results from coupling 
of configurations involving single-particle orbitals, {If, If) '■ ($fdf), {PfPf) and {dfdf), where 
the latter is already negligible at the present intensity. By selectively switching off final 
states consisting of (sfdf) and (PfPf) configurations we find that the IP ionization peak is 
dominated by [PfPf) orbitals while the CI ionization contribution is dominated by (s/o?/) 
contributions. These final states can be reached by absorption of two photons along the 
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LOPT pathways that correspond to either diagram (Fig. 2b) 



or diagram (Fig. 2a) 



xuv xuv , N , , 

Sf,Si — > p — > df) (20a) 



, xuv xuv , N f n , s 

(Pi -> s f ,pi -> df) (20b) 



While the first path (Eq. (|20ajl ) can be realized for the dominant configuration in the initial 
state (sj, Si) (see Eq. |£J)) and would be present for an uncorrelated initial state described by 
e.g. HF wavefunction, the second path (|20b|) has as prerequisite configuration admixtures 
(pi,Pi) to the initial state and thus initial-state angular correlation. When selectively elim- 
inating the (pi,Pi) configuration from the initial state we find that the cross section in the 
valley region is reduced by almost an order of magnitude. This unambiguously characterizes 
the "anomalous" cross section component in the "valley" as being due to correlations. By 
contrast, the IP ionization peaks are barely affected when (piPi) configurations are removed. 
This is plausible as the dominant two-photon absorption process from an uncorrelated initial 
state according to Fig. 2a 

, XUV XUV v ( . 

{Si -> p f ,Si -> pf) (20c) 

predicts a dominance of a Hertz dipole pattern for each ejected electron. Our calculations 
suggests that initial-state correlations may be more important than final-state correlations. 
This is due to the fact that the pair of electrons near the symmetric energy sharing point 
E s = 52 eV leave the interaction region quickly with a relatively large speed of v — 2 a.u. 

Ishikawa et al. Q| have discussed the "anomalous" component in terms of two semi- 
classical models. Post-ionization energy exchange (PIEE) and second ionization during core 
relaxation (SICR). They found that PIEE is inefficient to account for the valley region 
consistent with our observation that final-state correlations are of minor importance. On 
the other hand the relaxation process due to change in screening in the SICR appears to 
resemble somewhat a shake process and is as a quasi-isotropic process unlikely to yield a 
high-order (k = 4) anisotropy. 



V. CONCLUSIONS 



We have studied the electron energy and angular distributions in two-photon double- 
ionization of He by an attosecond, intense soft X-ray pulse, specifically, for the Thsapphire 

15 



59th harmonic pulse with an intensity of 10 15 W/cm 2 . We solved the TDSE with our coupled 
channel method in which the electron-electron interaction is fully taken into account. 

The electron energy distributions show well-localized peaks for pulse of long duration 
T p . They are understood to arise from the independent particle (IP) ionization. For short 
pulses of only a few hundred attoseconds, the peaks shift toward each other and the cross 
section in the valley between the peaks becomes significant. We attribute this ionization 
component to the correlation-induced (CI) ionization. We investigated the electron angular 
distributions from IP and CI ionization. We find shape profiles to be that of a Hertz dipole 
for IP ionization but a significant admixture of a k = 4 ("quadrupole") components for CI 
ionization. The unique signature of correlation-induced ionization is the presence of this 
k = 4 component in the angular shape profiles. They were further quantified in terms of the 
multipole expansion parameters. 

Time evolution of the electron angular distribution suggests that sequentiality of electron 
ejection or photon absorption is neither relevant nor well-defined. Clearly, further studies 
are needed to clarify electron correlation effects. Joint energy-angular distributions (i.e., 
kinematically complete momentum distributions) would provide new insight into the ion- 
ization mechanism. It would also be useful to understand the ionization with the help of 
a perturbative approach, either with the electron-electron interaction or the pulse inten- 
sity as the expansion parameter. This would provide a complementary picture to various 
mechanisms that may be difficult to identify in fully numerical TDSE results. 
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